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Abstract 

The theory for large amplitude circularly polarized waves propagating along an 
external magnetic field is extended in order to include also vacuum polarization effects. 
A general dispersion relation, which unites previous results, is derived. 

Circularly polarized waves can be exact solutions of the fluid and Maxwell equations 
(e.g. Refs. [1,2]). Thus, it is possible to study large amplitude waves propagating along an 
external magnetic field in a plasma, taking also relativistic effects into account. Here we 
shall demonstrate that it is comparatively easy to include quantum electrodynamical field 
(QED) effects [3, 4] in such an analysis, also for the general case with a multicomponent 
plasma. 

Let us start with the equations of continuity and momentum for each species (denoted 
by index s) in a cold, multi-component, relativistic, magnetized plasma 

d t N s + V • (N s v s ) = 0, (1) 

and 

d tPs + v s • V Ps = q s (E + v s x B) . (2) 

Here iV is the particle density, v the velocity, q the charge, m the rest mass, p = rav/(l — 
v 2 /c 2 ) 1 / 2 the momentum and c the speed of light in vacuum. Furthermore, we use the 
Maxwell's equations for the electric and magnetic fields E and B, i.e. 

V x E = -d t B, (3) 

and 

V X B - C~ 2 <9 f E = N sQs^s + W)Jvac, ( 4 ) 

s 

where the vacuum current J vac (after having corrected an obvious misprint in the definition 
of the coupling constant in Ref. [3]) is given by (e.g. Refs. [5,6]) 

J vac = -IkeI {c 2 V x [2{E 2 - c 2 B 2 )B - 7(E • B)E] 

-d t [2(£ 2 -c 2 £ 2 )E + 7c 2 (E-B)B]}, (5) 

where k = 2a 2 h 3 /45mfc 5 is the nonlinear coupling constant and a = e 2 /^ttTiceq is the fine 
structure constant. Here e (= — q e ) is the elementary charge, m e the electron rest mass, 
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£q the vacuum permittivity, //o the vacuum permeability and h is 2tt times the Planck 
constant. The expression J5)) for the vacuum current is valid for field strengths lower than 
the critical field strength, E cr i t = m 2 c s /he, and for frequencies lower that the Compton 
frequency, uj e = m e c 2 /h. These equations have been recently used by Di Piazza et al. [6]. 
We now consider the propagation of a circularly polarized wave along an external 



constant magnetic field Bq z z. 
then constants and that E • B 



Following Ref. [2] we first note that E as well as B are 



0. This means that 



J vac = -Ake 2 q (E 2 - c 2 B 2 ) (c 2 V x B - dtE) . 
Inserting Eq. J6]) into Eq. (|4|) we thus find our new equation 

V X B - C~ 2 dtE = /i Jeff, 



where 
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Here n = kc/u is the index of refraction, and E ± is the electric field amplitude of our 
right hand circularly polarized wave E = Eq± [cos (cot — kz) x + sin [uit — kz) y] . The rest 
of the calculations are now straightforward. Introducing the symbols E± = E x ± iE, 
B± = B x ± iB y and v ± = v x ± iv y we thus adopt the analysis of Ref. [2], noting that 

E± = E ± exp (±iuit =f ikz) , 



B± = ±i-E±, 



(8b) 



and 



7 s m s (oj + uj cs ) 



E±, 



where uj cs = qsBo z /j s m s is the relativistic gyrofrequency of particle species s and 7 = 
(1 — Vq/(?)~ 1 / 2 is the Lorentz factor which in this case is a constant with Vq = v+v_. 
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Usually, we define the relativistic plasma frequency as uj ps = (q 2 No s /~/ s m s eo) . Here, 
however, instead of uj ps we will in the dispersion relation find the effective plasma frequency 
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By combining Eqs. (|7| and |(SJ) we obtain the dispersion relation 
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(n 2 - 1) . (9) 



The theory for a left hand circularly polarized wave is quite similar and leads to the same 
dispersion relation but with oj cs — ► — u> cs . 
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We note that Eq. ([9]) agrees with Eq. (13) in Ref. [7] in the special case of an ultra- 
relativistic electron-positron plasma, and with Eq. (9) in Ref. [8] for a dusty plasma (if we 
correct for a misprint in Refs. [7] and [8]; 4a should be 2a). Eq. j9]) is exact within the 
basic model we have used in the present paper. 

As aE 2 ± / E 2 rit in general is a very small parameter, we can next expand the solution 
of Eq. © to obtain 
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(10) 



which agrees with the result of Ref. [6] in the limit Bq z = 0. 

New modes in unmagnetized QED plasmas have been reported previously [9]. The 
calculation were then extended to magnetized QED plasmas with ultrarelativistic ions. 
Such plasmas can support amplitude-dependent electromagnetic ion waves [10] 
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where u> p i is 



the ion plasma frequency and ooEi = e£oi/cwij, as well as new modes in 
a dust-ion plasma [11], and QED effects can then play a key role [8, 11]. Di Piazza et 
al. demonstrated in Ref. [6] that the QED effects can be important in wave propagation 
studies, even for the case where Bq z = 0. They showed that the last term in Eq. (TlOl) 
can be larger than the sum of the first two terms on the right hand side of Eq. (jTOl) 
near the plasma cut-off frequency. Thus, they demonstrated enhancement of the vacuum 
polarization effects when the frequency is close to lo p . Our dispersion relation shows that 
there is no analogous enhancement of the QED contribution from the external magnetic 
field near the plasma cut-off frequency. This is because the 2nd and the 3rd term on the 
right hand side of Eq. (fTOj) have the same sign. Thus, the QED contribution from the 
external magnetic field will always remain smaller than the classical contribution in Eq. 
(fTOl) unless cBq z 3> E CT n, but this condition lies outside the validity of our model. The 
QED contribution from the external magnetic field may still be larger than that of the 
last term in Eq. (fTOl) . however, in particular if cBq z > Eq±, lo p < lo and if (lo + to cs ) is not 
close to resonance. But in this regime there is no enhancement of the vacuum polarization 
effects as described in Ref. [6]. 

It is also interesting to note that the effective plasma frequency Q p is different from to p . 
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However, for field strengths available in laboratory experiments, this ratio is close to unity. 

Another example of interest concerns the free electron laser (FEL) [12-14]. Considering 
magnetostatic wigglers, i.e. letting lo — > 0, in a one-component electron plasma with drift 
velocity v ze and gyrofrequency to ce , the dispersion relation using the laboratory frame 
rather than the electron rest frame, becomes 



where we have introduced the total magnetic field B tot = y B^ ± + B^ z and the magnetic 
field amplitude B j_ of the wiggler field. For cB to t *C E CI n, Eq. (fTTj) agrees for example 
with (8) in Ref. [14]. Thus our formula (fTTj) generalizes the FEL dispersion relation to ac- 
count for field strengths approaching the Schwinger critical field. While current laboratory 
applications have a long way to go in order to have fields in this range, there nevertheless 
exist schemes attempting to reach such regimes [15]. We note that the particle orbits in 
FELs are nontrivial [13,14] and that the helical equilibria may not be stable [13]. In the 
example above we have also supposed the presence of a stationary background assuring 
neutrality, which means that there are several limiting assumptions concerning forthcom- 
ing experiments [6]. Thus we hope that the readers will find interest in further extending 
the basic theory of Ref. [2]. 

In the present manuscript we have derived the dispersion relation for large amplitude 
circularly polarized electromagnetic waves propagating parallel to an external magnetic 
field in a multicomponent plasma. Our dispersion relation (j^J) accounts for vacuum po- 
larization and magnetization (both due to the external field and the wave field), fully 
relativistic quiver velocities and general drift velocities of the plasma species, and applies 
for arbitrary ratios of the characteristic frequencies of the problem. As a consequence, 
such different waves as Alfven modes, whistler modes and large amplitude laser modes are 
covered by the general theory. Considering specific regimes, Eq. (JjJJ) reduces to several 
special cases found in the literature. Besides uniting the results of many previous works in 
a single formalism, the results presented here can find applications to wave propagation in 
astrophysical contexts, e.g. in magnetar atmospheres [16], as well as in the next generation 
of FELs [15]. 

We consider the present work as a necessary prerequisite to a more complete theory 
which also includes a full stability analysis [2]. Thus it has for example been shown [1] 
that a single cold electron beam in an unmagnetized plasma is unstable if k = and if 



where and K are the frequency and wavenumber of the perturbations. The vacuum 
current will however modify all stability criteria. 
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